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ABSTRACT 



We present Herschel SPIRE observations for the TW Hydrae association (TWA) brown dwarf discs SSSPM 11102-3431 (SSI 102) 
and 2MASSW J 1207334-393254 (2M1207). Both discs are undetected in the SPIRE 200-500yi/m bands. We have also analyzed the 
archival PACS data and find no detection for either source in the 160/jm band. Based on radiative transfer modeling, we estimate an 
upper limit to the disc mass for both sources of 0.1 M ]up . The lack of detection in the SPIRE bands could be due to a paucity of 
millimeter sized dust grains in the 2M1207 and SSI 102 discs. We also report a non-detection for the brown dwarf 2MASS J 1 1 395 11- 
315921 (2M1 139) in the PACS 70 and 160/jm bands. We have argued for the presence of a warm debris disc around 2M1 139, based 
on an excess emission observed at 24//m. The mid-infrared colors for 2M1139 are similar to the transition discs in the Taurus and 
Ophuichus regions. A comparison of the brown dwarf disc masses over a ~1-10 Myr age interval suggests a decline in the disc mass 
with the age of the system. 
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y 1. Introduction 

^2 ! The TW Hydrae association (TWA) has been under observation 
i i ■ since its identification by Kastner et al. (1997). This young (10*5 
Myr; Barrado y Navascues 2006) and nearby (~53pc, Gizis et al. 
2007) association now has more than 30 confirmed members, in 
both stellar and sub-stellar regimes, with a few members dis- 
playing signatures of being in the T Tauri phase. Given its prox- 
imity and a critical intermediate age when most discs begin to 
show a transition to the debris phase, TWA is well-suited for 
the study of disc evolution and planet formation. Using ground- 
based mid-infrared observations, Jayawardhana et al. (1999) and 
Weinberger et al. (2004) were able to detect excess thermal emis- 
sion at 10, 12 and 18 pm for four TWA members (TW Hya, Hen 
3-600, HD 98800A and HR 4796A). However, most of the other 
members do not show any evidence for circumstellar dust at 
these wavelengths, indicating cleansing of dust in the inner disc 
regions on a time scale of <10 Myr. A Spitzer survey by Low et 
al. (2005) using the sensitive MIPS instrument revealed the pres- 
ence of cold debris discs around two TWA members, TWA 7 and 
1 3, at 70 pm, that remained undetected by IRAS and other previ- 
ous surveys. Among the brown dwarf members of TWA, excess 
emission has been detected at mid- and far-infrared wavelengths 
for 2MASSW J1207334-393254 (2M1207) and SSSPM J1102- 
3431 (SSI 102), while a small 24 pm excess has been reported 
for 2MASSW Jl 13951 1-315921 (2M1 139), suggesting the pres- 
ence of a transition disc around this source (Riaz et al. 2006; 
Riaz & Gizis 2008). By transition discs, we imply a disc with an 
inner opacity hole. The inner disc in such cases may have opti- 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



stars: low-mass, brown dwarfs - stars: individual (2MASSW J 1207334- 



cally thin dust, or could be completely devoid of circumstellar 
material. Due to reduced opacity in the inner disc regions, tran- 
sition discs exhibit photospheric emission shortward of a certain 
mid-infrared wavelength, depending on the size of the inner disc 
hole. Based on the overall disc fractions and the strength in the 
mid-infrared excess emission, Riaz & Gizis (2008) suggested 
longer disc dissipation timescales for the brown dwarfs in TWA. 
This association has proven to be an important field to study the 
formation and evolutionary differences, if any, in stellar and sub- 
stellar objects, and their disc evolution time scales. 

We present ESA Herschel Space Observatory (Pilbratt et al. 
2010) SPIRE (Griffin et al. 2010) observations at 200-500/an for 
SSI 102 and 2M1207, and PACS (Poglitsch et al. 2010) 70 and 
160/im observations for 2M1139. Observations at far-infrared 
wavelengths for SSI 102 and 2M1207 can help constrain the 
model fits and obtain better estimates for the disc mass and 
outer disc radius, while PACS observations of 2M1139 are im- 
portant to probe the presence of cold dust material for this tran- 
sition disc. Our Herschel observations are presented in Section 
§|2] Results from the radiative transfer modeling are presented in 
Section §|3] Section $4]provides a comparison of the TWA brown 
dwarf disc masses with the higher mass stars in this association 
and with the younger brown dwarf discs in Taurus, along with a 
discussion on the nature of the transition disc around 2M1 139. 



2. Observations and Data Analysis 

We obtained SPIRE photometric observations for SSI 102 and 
2M1207 in the 250, 350 and 500/i/m bands, using the small map 
mode (ObsID: 1342234825/1342223257). The small map mode 
covers a 5' diameter circle, with a fixed scan angle of 42.4° and 
a fixed scan speed of 30"s~' . We requested a setup of 25 repeti- 
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tions, resulting in a total on-source integration time of 925s. For 
2M1139, PACS photometric observations were obtained in the 
scan map mode in the blue (10pm) and the red (160/im) filters. 
Two different scan map angles of 45° and 63° were used (ObsID: 
1342224188/89). Both scans were obtained at a medium scan 
speed of 20"/s, with a cross scan step of 5" and scan leg length 
of 3'. The number of scan legs was 8 for the 63° scan, and 4 for 
the 45° scan. 

Data analysis for all sources was performed using the rectan- 
gular sky photometry task provided in the Herschel Interactive 
Processing Environment (HIPE version 8.0.1). We have worked 
on the final pipeline processed data (the 'Level 2' products). 
The rectangular photometry task applies aperture photometry 
for a circular target aperture and a rectangular sky aperture. 
For 2M1139, photometry was performed on the mosaic of the 
two scan maps. The sky intensity was estimated using the me- 
dian sky estimation algorithm. Since the sky intensity could 
vary with location, we selected 4 different background regions 
around the target (each of 30x30 pixel area) and measured the 
sky-subtracted source flux for each sky region. The final source 
flux and the flux error is the mean and the standard deviation of 
these 4 flux measurements. The flux errors are consistent with 
the values obtained at the source location from the 'error' image 
provided in the Level 2 products. 

2.1. SS1102 

SSI 102 is undetected in all SPIRE bands (Fig. [TJ). The posi- 
tional accuracy as determined from the actual ra, dec and the ra- 
Nominal, decNominal in the metadata is ~ 1 pixel (or ~6" in the 
250/mi band). We have calculated the upper limits using a pixel 
radius of 2 pixels centered at the nominal position of SSI 102. 
Harvey et al. (2012) have reported a detection for SSI 102 in the 
PACS 160/mi band, at a ~7 mJy level. We have analyzed their 
PACS data, and find no source detection at the nominal or the 
actual position of the target (Fig. [2]>. The flux level at the target 
location is ~0. 1-0.2 mJy, and is the same as found in some of the 
confusion noise dominated regions seen in the image. The off- 
set between the actual and the nominal positions in the metadata 
is ~1.4", or 3.5 pixels. However, there is no source detection at 
a 2-cr level even if the offset is taken into account. The 160yum 
flux measurement for SS 1 102 reported by Harvey et al. (2012) is 
thus incorrect. Our analysis of SSI 102 in the PACS 70yum band, 
however, is consistent with the measurement reported by Harvey 
et al. (2012). Tableffllists the updated fluxes for SSI 102. 

2.2. 2M1207 

We had reported in Riaz et al. (2012a) a detection for 2M1207 
in the SPIRE bands of 250 and 350/vm, but a comparison now 
with the PACS images indicates this to be a misidentification. 
The revised SPIRE fluxes (upper limits) for 2M1207 have been 
presented in Riaz et al. (2012b). As noted in that work, there is a 
possibility of having an unresolved source in the SPIRE bands, 
given the large positional offset of ~14". However, the contam- 
ination from the nearby bright object is likely to dominate the 
composite photometry in such a scenario. The SPIRE upper lim- 
its have been determined by measuring the emission at the nom- 
inal position for 2M1207, without considering the offset. 

We have also analyzed the PACS 70 and 160/mi observations 
for 2M1207. Our 70/mi measurement is consistent with the flux 
reported by Harvey et al. (2012). In the 160yum image, there is 
some diffuse emission seen at a distance of 2 pixels (~0.8") from 



the nominal position of 2M1207 (Fig. [2]). This diffuse emission is 
at a flux level of ~0.2mJy, similar to that observed at the location 
of SSI 102 in the 160//m data. We can consider this flux level 
as the expected confusion noise in the 160/im band, based on 
which the detection of both 2M1207 and SSI 102 at 160jum can 
be rejected. 

2.3. 2M1139 

2M1139 is undetected in the PACS bands of 70 and 160jt/m 
(Fig-EJ- The position uncertainty is ~0.5 pixel (~0.2"). We have 
determined the upper limits by measuring the flux at the nomi- 
nal position of the target. We note that the Spitzer IQpm upper 
limit reported in Riaz et al. (2009) is much higher than the PACS 
10pm measurement, which can be expected given the higher 
confusion noise in the Spitzer MIPS bands compared to PACS. 



3. Disk Modeling 

Figure H shows the model fits for SSI 102 and 2M1207. The fits 
are based on 2-D radiative transfer modeling of near-infrared to 
sub-millimeter data, using the code by Whitney et al. (2003). 
Detailed discussion on the fitting procedure and the variations in 
the model SEDs is provided in Riaz & Gizis (2007). Here we 
have provided a brief description of the best model-fits obtained, 
based on the lowest x 1 value. There are mainly six parameters 
related to the disc emission which can be varied to obtain a good 
fit (Table 0. The disc scale height varies as h — ho{m / 'R*f ', 
where ho is the scale height at R„ and /3 is the flaring power. 
Both 2M1207 and SSI 102 show flat disc structures, and val- 
ues for /? of ~ 1.1-1.2, and ho of ~0.3 provide a good fit. We 
have used the large grain model (a m „=lmm) in the disc mid- 
plane and the medium sized grains (a max ~ 1 pm) in the disc 
atmosphere. Since both discs show no silicate emission at \0pm 
(Riaz & Gizis 2008), grains larger than the ISM are required to 
fit the flat silicate feature. Further details on the variations in the 
model SEDs for the different grain models are discussed in Riaz 
& Gizis (2007). Due to binning of photons in the models, there 
are a total of 10 viewing angles available. For 2M1207, a good fit 
was obtained using close to edge-on inclination angles between 
57° and 69°, whereas SSI 102 is modeled well using intermedi- 
ate angles of 53°-63°. We have varied the inner disc radius R m i„ 
in multiples of R su b, which is the dust sublimation radius. We 
have obtained good fits for R m j n between 1 and 3R su i„ with a 
peak at lR su b (based on a^- 2 value comparison). This implies an 
absence of an inner hole in the disc since it would have to be 
larger than the dust sublimation radius. Using the sub-mm upper 
limits, the best model fit for both 2M1207 and SSI 102 is for an 
outer disc radius, R max , of 50 AU and an upper limit to the disc 
mass of 0.1 Mj up . A 100 AU model for a mass of 0.5 Mj up is 
also a good fit, though at a slightly higher^ 2 value than the best 
fit. The 10 and 30 AU models are not a good fit to the longer 
part of the IRS spectrum for wavelengths >l0pm. The optical 
depth is too large for these models even if the disc mass is re- 
duced to 0.01 Mj up , and misses the 70-350/mi points. Increasing 
the disc mass for the same radius increases the optical depth be- 
yond ~30yum, while increasing the radius for the same disc mass 
reduces the optical depth longward of ~30pm. The best fits for 
2M1207 and SSI 102 are the same as provided in Riaz & Gizis 
(2008). We refer the readers to that paper for details on the rest 
of the fitting parameters. 

We have also attempted to model the 24pm excess along 
with the 70 and 160/mi upper limits for 2M1139 (Fig. g). We 
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Fig. 1. The SPIRE images for SSI 102. Blue circle marks the actual target position, green circle marks the nominal position. The 
image scale for 30" is shown in the top right. North is up, East is to the left. 
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Fig. 2. The PACS 160//m images for SSI 102 (left) and 2M1207 (right). Blue circle marks the actual target position, green circle 
marks the nominal position. For the SSI 102 image, black circles indicate some of the confusion noise dominated regions with 
similar flux levels as observed at the nominal target position. The image scale for 10" is shown in the top right. North is up, East is 
to the left. 



have used the large grain model (a max =\mm) in both the disc 
midplane and the atmosphere to fit this source. The best fit is 
obtained for a disc mass of 1E-7 Mq and an outer disc radius 
of 100 AU. Since this is a transition disc with an inner opacity 
hole, a larger value of R m i„~l R su b provides a good fit. The pa- 
rameter values listed in Table |2]for 2M1 139 model fit should be 
considered with caution. These or any of the other disc param- 
eters cannot be properly constrained, considering that we have 
only one excess point and two upper limits. This is particularly 



for the case of the disc inclination angle and the outer disc ra- 
dius. The model fit for 2M1 139 mainly indicates that there may 
be a small mass of warm dust located at ~ 1 AU distance from 
the central brown dwarf. 

4. Discussion 

Figure [5] shows the extent of excess emission at 24 /mi for the 
TWA stars and brown dwarfs. The dot-dashed line represents 
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Fig. 3. PACS images for 2M1 139. Blue circle marks the actual target position, green circle marks the nominal position. The image 
scale for 10" is shown in the top right. North is up, East is to the left. 



SSSPM 1102 




-0.5 0.5 1 1.5 2 2.5 



log X (cm) 




log X (urn) 

2M1139 




log X (urn) 



Fig. 4. The best model fits for SSI 102, 2M1207 and 2M1139 (bottom). Also shown is the contribution from the disc (blue) and the 
stellar photosphere (grey). The Spitzer/IRS spectrum is shown in red. Blue points are the Herschel upper limits. For 2M1 139, the 
red point marks the Spitzer 70pm upper limit. 



the limit of FiaIFk under the assumption that both bands lie 
on the Rayleigh- Jeans tail of the stellar spectrum. The F 24-to- 
F K ratio for 2M1139 is higher than TWA 7 or TWA 13, indi- 
cating a warmer debris disc around this brown dwarf compared 
to the cold 80 K disc around TWA 7 or the 65 K disc around 
TWA 13AB. The dashed line represents a geometrically thin, op- 
tically thick flat disc with a spectral slope AF^ cc A~ 4 ^ . 2M1207 
and SSI 102 lie just at the dashed line, indicating optically thick 
flat discs around these two brown dwarfs. The primordial na- 
ture of 2M1207 and SSI 102 discs has also been argued in Riaz 
et al. (2012c), based on the strength in the excess emission ob- 



served in the Spitzer IRAC bands and the IRS spectrum. That 
these discs are optically thick can also be inferred from the K s 
to 24-pm slope, a, of the SED. We can consider a boundary of 
a = -2.2 to separate the Class II and Class III sources (e.g., 
Luhman et al. 2010). 2M1207 and SSI 102 have a of -1.38 and 
-1.42, respectively, consistent with being Class II systems, while 
2M1 139 has a of -2.17 and is just above the Class III boundary. 
Large negative spectral slope values can be expected for transi- 
tion discs with photospheric emission shortward of 24/im due to 
inner opacity holes (e.g., Cieza et al. 2010). Transition discs can 
be better characterized in a a excess vs. A turno ff diagram (Pig. [6j, 
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where Atumoff is the longest wavelength at which the observed 
emission is photospheric, and a excess is the spectral slope be- 
tween A lurno ff and 24//m (e.g., Cieza et al. 2010). 2M1139 has 
Atumoff of 8jum and a excess of -2.1. The parameter A tunm ff is a 
good indicator of the size of the inner opacity hole in a tran- 
sition disc. For a sub-stellar source, A, urno ff of 8jL/m implies an 
inner hole of ~1 AU (e.g., Muzerolle et al. 2006; Cieza et al. 
2010). We have included in Fig. [6] the sample of young transi- 
tion discs in the Ophichus and Taurus regions presented in Cieza 
et al. (2010; 2012). A large negative value for a excess indicates 
a debris disc or a continuous disc that has undergone significant 
grain growth and dust settling, compared to large positive val- 
ues that indicate an outer optically thick disc with a sharp inner 
cut due to e.g., the presence of a companion in the inner opac- 
ity hole. The location of 2M1139 in Fig. [6] is consistent with 
being a warm debris disc. It has a low fractional disc luminosity 
(Ld/L„<10~ 4 ) and is a non-accretor(< 10 _11 MQyr _1 ), similar to 
the debris disc candidates in the Ophichus and Taurus regions. 

Fig.|7]compares the mid-IR colors for the TWA brown dwarf 
discs with the Class II and the discless sources in Taurus (from 
Luhman et al. 2010) and the transition disc samples in Ophichus 
and Taurus from Cieza et al. (2010; 2012). The fluxes at 3.6 
and 4.5/mi are not known for SSI 102 and so it cannot be in- 
cluded in the top panel in Fig. [7] The dashed lines are boundaries 
used by Cieza et al. (2010) to separate the transition discs from 
the "full" discs with no inner holes. The colors for 2M1 139 are 
similar to the sources classified as debris discs in these regions 
(black points), while 2M1207 and SSI 102 lie among the Class 
II sources. Transition discs with flattened disc structures due to 
grain growth/dust settling show redder [3.6]-[24] and [5.8]-[8] 
colors than the debris discs. While 2M1207 and SSI 102 lie to 
the right of the transition disc boundary, their colors are also 
similar to discs dominated by grain growth/dust settling. This is 
consistent with the flat 10 and 20jt/m silicate emission features 
observed for both 2M1207 and SSI 102 (Riaz & Gizis 2008), 
indicating grain growth to sizes >2jum followed by dust set- 
tling to the disc midplane. We can expect weaker excess in the 
[3.6]-[24] or [5.8]-[8] color for a -10 Myr old disc compared to 
the younger ~1 Myr Class II sources in Taurus, most of which 
are still dominated by sub-micron sized dust grains (Riaz 2009; 
Sargent et al. 2009). 

There are four widely discussed mechanisms to explain the 
formation of inner opacity holes in transition discs. The process 
of giant planet formation would require a large disc mass (> 1 
Mj up ) and therefore can be rejected for 2M1139. Tidal trunca- 
tion due to a close companion may also not be applicable for 
2M1139, as Chauvin et al. (2010) have found no secondaries 
for this object at projected separations of ~5-240 AU. However, 
given an inner hole size of ~1 AU for 2M1 139, we cannot rule 
out the possibility of a tight binary surrounded by a circumbi- 
nary debris disc. We note that there is no report that 2M1139 
is a close double in the limited set of published high resolution 
spectra. Inner disc evacuation due to significant grain growth and 
dust settling is expected to exacerbate accretion, resulting in an 
actively accreting disc with an inner opacity hole (e.g., Cieza 
et al. 2010). This mechanism would also not be applicable for 
2M1 139. The transition from a primordial to a debris stage could 
take place due to the photoevaporation mechanism, which can 
set in once the accretion rate drops below the photoevaporative 
wind rate. As soon as photoevaporation sets in, the remaining 
gas in the disc is removed very rapidly (<< 1 Myr), and what 
remains is a gas-poor debris disc (e.g., Alexander & Armitage 
2006). For discs that are still in the process of photoevaporation, 
we can expect weak or negligible accretion but larger fractional 



disc luminosities than 2M1 139, of the order of 10~ 2 (Fig. [6]). The 
much lower L^/L, for 2M1 139 thus suggests that it has already 
made a transition to the debris phase. Interestingly, while the 
24/vm excess for 2M1 139 is larger than that observed for the de- 
bris discs TWA 7 and 13 (Pig. [5j, it remains undetected at 70/im 
and beyond. This is unlike TWA 7 and 13 which show strong 
excesses at 70/mi (Low et al. 2005). This suggests an absence 
of cold dust material for 2M1 139, and perhaps the outer disc re- 
gion has been dissipated completely. We cannot, however, con- 
clude that there is no cold dust material surrounding 2M1139, 
or that this is a case of inner disc clearing due to photoevapora- 
tion. There is a limit to the fractional disc luminosity that can be 
reached with a certain exposure time. We estimate that we can 
detect an AU Mic like disc (L D /L„ ~ 6E-4; Liu et al. 2004) with 
the ~1 hour total exposure time used for 2M1 139. The fractional 
disc luminosity as estimated from the 2M1 139 SED is of the or- 
der of ~ 10~ 5 , which is well below the detection limit that can 
be obtained with the requested setup. 

There are results suggesting two distinct populations in the 
TWA, as evidenced by a bimodal distribution in the rotation pe- 
riods (Lawson & Crause 2005) and the presence of warm (T^ 
100 K) dust (Low et al. 2005; Weinberger et al. 2004). Low 
et al. found negligible amounts of warm dust around 20 out of 
their 24 TWA targets, while the other four display strong excess 
emission at 24 //m (Fig. [5]). Lawson & Crause have found the 
projected rotational velocities, v sin i, for the TWA 1-13 group 
led by TW Hya to be 3-10 km s , while the velocities for the 
stars in the TWA 14-19 group are ~20-50 km s _1 . These authors 
have suggested that the former group might be younger than the 
latter by ~8 Myr (8-10 Myr versus -17 Myr). 2M1139 has a v 
sin i of 25 kirT 1 , which is a factor of ~2 higher than 2M1207 
(13 km s _I ; Scholz et al. 2005). This brown dwarf thus could be 
slightly older than 2M1207, which could explain its advanced 
evolutionary stage of being a debris disc, compared to the pri- 
mordial colors observed for 2M1207 and SSI 102. We note that 
a variety in SED morphologies is observed at any given age (e.g., 
Hartmann et al. 2005), and the presence of a debris disc does not 
necessarily imply an older age. 

2M1207 is known to have a ~5Mj up secondary that lies at a 
projected separation of ~55 AU (Chauvin et al. 2004). The best 
model fit for 2M1207 is obtained for an outer disc radius of 50 
AU, which suggests tidal truncation of the disc by the massive 
companion. Such a tidal truncation could be responsible for the 
sub-mm non-detection of this disc. However, we cannot rule out 
the non-detection for both 2M1207 and SSI 102 in the SPIRE 
bands due to limitations by the confusion noise. We have dis- 
cussed in Riaz et al. (2012a) the possibility of the formation of 
the ~5Mj up secondary in the 2M1207 disc via disc fragmenta- 
tion. Our argument is that even a very low mass disc could pro- 
duce a fragment of ~Q.Q35Mj up , which can then grow over time 
to form a 5 Mj up mass object. The main requirement for such a 
case is for the initial mass of the disc to be higher than its cur- 
rent estimate (at least 10-20 Mj up ). An upper limit on the disc 
mass of ~0.1 Mj up as obtained from modeling of 2M1207 SED 
thus does not rule out disc fragmentation, since the system is 
relatively old and we have not observed it during its early stages 
when fragmentation could have occurred (<0.1 Myr). This sys- 
tem can rather be considered as a wide very low-mass binary 
than a planetary system, and some alternative formation mecha- 
nisms, such as turbulent fragmentation of a core, have been dis- 
cussed in Riaz et al. (2012a). 

Figure [8^ compares the relative disc mass upper limits for 
2M 1207, SSI 102 and 2M1 139 with the other disc sources in the 
TWA (using data from Calvet et al. 2002, Low et al. 2005 and 
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Fig. 5. The extent of excess emission at 24 /urn for the TWA stars and brown dwarfs. The dot-dashed line represents the limit of 
Fia/Fk under the assumption that both bands lie on the Rayleigh- Jeans tail of the stellar spectrum. The dashed line represents a 
geometrically thin, optically thick flat disc with a spectral slope AF A oc A~ 4/i . 
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Fig. 6. A a excess vs. A, umo ff diagram showing 2M1139 and the transition disc samples in Ophuichus and Taurus. Transition discs 
with inner holes due to grain growth/dust settling are indicated by red symbol, due to photoevaporation by blue, and due to giant 
planet formation by grey. Debris discs are indicated by black filled symbol. Open green circle indicates objects with accretion rates 
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Lq/L, between 1E-3 to 1E-2.. Black open circles indicate non-accretiors (<1E-1 1 Mq yr~') with weak discs (Lo/L* < 1E-3). 
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Fig. 7. Mid-IR color-color diagram for TWA brown dwarf discs. Also included for comparison are the discless (grey) and Class II 
(blue) sources in Taurus, and the transition disc sample in Ophuichus and Taurus. Green symbol indicates the transition discs with 
inner disc evacuation due to grain growth/dust settling. Black symbols indicate the debris discs. 



Andrews et al. 2010). TW Hya is known to be a strong accretor, 
with M of 4 x 10" 10 M o yr l (Muzerolle et al. 2000). Hen 3- 
600 is a weaker accreting source (0.5xl0" 10 M© yr- 1 ), while the 
other disc sources in TWA are non-accretors (e.g., Muzerolle 
et al. 2000). Among the brown dwarfs, 2M1207 shows vari- 
able accretion activity of 0.8-1.6 x 1O~ IO M0 yr' 1 (Scholz et 
al. 2005; Stelzer et al. 2007), while SSI 102 and 2M1139 are 
non-accretors (< 10~" Mq yr' 1 ; Scholz et al. 2005). The disc 



mass limits for 2M1207 and SSI 102 are comparable to the 
weak/non-accretors in the TWA, while 2M1 139 is more similar 
to the other two transition discs TWA 7 and 13. Transition discs 
are found to have much lower disc masses and accretion rate 
than non-transition disc sources (e.g., Cieza et al. 2010). Both 
accreting/non-accreting objects, or sources with or without mid- 
IR excess are found to have an equal likelihood of a sub-mm de- 
tection (e.g., Lee et al. 201 1). The lack of detection in the SPIRE 
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bands for 2M1207 and SSI 102 could be due to a paucity of mil- 
limeter sized dust grains in these discs. Fig.lHJ) compares the rel- 
ative disc masses for the Taurus, IC348 and TWA brown dwarfs, 
thus covering an age range of ~1-10 Myr. The Taurus and IC 
348 disc mass measurements are from Scholz et al. (2006) and 
Klein et al. (2003), respectively. The TWA discs have relative 
masses smaller than nearly all brown dwarf discs in Taurus, thus 
suggesting that the disc mass may be dependent on the age of 
the system. For higher mass stars (~0.1-1 Mq) in the Taurus, 
Ophiuchus and IC 348 clusters, the disc mass distribution for 
the relatively older 2-3 Myr clusters is found to be shifted by 
a factor of ~20 to lower masses, compared to the distribution 
for the ~1 Myr old clusters (Lee et al. 2011). For the case of 
the brown dwarfs, most of the data points in Fig. [8j> for object 
masses < 0.03 Mq are upper limits. The non-detection of these 
faint objects in the sub-mm/mm bands could be due to low sen- 
sitivity of the observations. To probe any age dependence of the 
disc mass would require larger samples of brown dwarfs across 
a wide range in ages. 

5. Summary 

We report a non-detection of SSI 102 and 2M1207 in the SPIRE 
and PACS 160yum bands, and a non-detection of 2M1 139 in the 
PACS 70 and 160//m bands. 2M1139 likely possesses a warm 
debris disc, as indicated by its 24/im excess emission. We esti- 
mate an upper limit to the disc mass for SSI 102 and 2M1207 
of 0. 1 Mj up , which places them among the weaker brown dwarf 
discs in the ~1 Myr old Taurus region, and suggests a decline in 
the brown dwarf disc mass with the age of the system. 
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Fig. 8. Top- (a): Relative disc masses vs. object mass for TWA discs. Bottom: (b) A comparison of TWA relative disc masses with 
other known brown dwarf discs in Taurus and IC 348. Upper limits are shown as bold triangles. 
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Table 1. Observations 



Band" 


2M1207 [Ref] 


SSI 102^ [Ref] 


2M1139 [Ref] 


/ 


1.13+0.1 [1] 


0.83+0.1 [1] 


1.17+0.1 [1] 


./ 


10.10±0.89 [2] 


9.75+1.0 [2] 


13.43±0.3 [2] 


H 


11.35±1.0 [2] 


1 1.7+1.0 [2] 


16.3±0.3 [2] 


K 


11.12±0.98 [2] 


11.72+1.0 [2] 


16.7±0.3 [2] 


3.6yum 


8.49±0.32 [3] 


— 


12.08+0.7 [3] 


4.5yum 


7.15±0.26 [3] 


— 


8.5±0.6 [3] 


5.8 fim 


6.36±0.06 [3] 


7.6+0.2 [6] 


5.98±0.6 [3] 


8yum 


5.74±0.21 [3] 


6.74±0.2 [6] 


3.88+0.5 [3] 


3.8yum 


7.00+0.6 [5] 


— 


10.23±1 [5] 


8.7yum 


3.0U+1 [Oj 






10.4yum 


7.50±1 [6] 






24yum 


4.32±0.03 [3] 


3.96±0.1 [3] 


1.10+0.2 [3] 


70/jm 


9±4 [4] 


9.4±0.6 [7] 


<0.2 [8] 


160/im 


<0.2 [8] 


<0.2 [8] 


<0.075 [8] 


250/jm 


<5.2 [8] 


<1.5 [8] 




350yum 


<5 [8] 


<1 [8] 




500yum 


<16 [8] 


<1 [8] 





a Fluxes in [mjy]. References: 1 - DENIS; 2 - 2MASS; 3 - Riaz et al. (2006); 
4 - Riaz & Gizis (2008); 5 - Jayawardhana et. al. (2003); 6 - Sterzik et. al. 
(2004); 7 - Harvey et al. (2012); 8 - This work. 

* The fluxes at 5.8, 8 and 24/um are from the Spitzer/IRS spectrum. 



Table 2. Disk parameters 


Parameter 


2M1207 


SSI 102 


2M1139 




1.2+0.01 


1.12+0.01 


1.0+0.01 


h) 


0.3±0.1 


0.25+0.1 


0.1+0.1 




1 Rsub (~3 Rm) 


1 Rsub 


7 Rsub 




50 - 100 AU 


50 - 100 AU 


100 AU 


M disc 


<0.1 M Jup 


<0.1M,„ P 


<lE-7 M Jup 


i 


57° - 69° 


53° -63° 


50° 



